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We predict through our first-principles calculations that four double perovskite oxides of Bi2ABC>6 (AB 
= FeMo, MnMo, MnOs, CrOs) are half-metallic ferrimagnets. Our calculated results shows that the four 
optimized structures have negative formation energy, from -0.42 to -0.26 eV per formula unit, which implies 
that they could probably be realized. In the case of Bi 2 FeMo06, the half-metallic gap and Curie temperature 
are predicted to reach to 0.71 eV and 650 K, respectively, which indicates that high spin polarization could be 
kept at high temperatures far beyond room temperature. It is believed that some of them could be synthesized 
soon and would prove useful for spintronic applications. 

PACS numbers: 75.30.-m, 75.50.-y, 71.20.-b, 75.10.-b 
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I. INTRODUCTION 

Magnetic materials that have high spin polarization at 
room temperature or higher are highly desirable for spin- 
tronic applications 1 ' 2 . Half- metallic materials are good 
candidates because they can have high spin polarization 
at high temperature 3 . In the case of CrC>2, 96% spin 
polarization has been achieved experimentally 2 . Since 
1998, double perovskite oxides have been explored ex- 
tensively for this purpose because both half-metallicity 
and high Curie temperature can be achieved in such 
materials^. It has been shown experimentally that sev- 
eral double perovskite oxides, such as S^FeMoOe and 
S^CrReOg materials^, can keep their ferromagnetic or 
ferrimagnetic phases far beyond room temperature, and 
more importantly, high-quality materials have been real- 
ized recently^— . 

Here, we present our first-principles exploration on 
double perovskite oxides of Bi2ABC>6 with A being 3d 
transition metals and B 4d/5d ones. We optimize their 
structures fully and then investigate their stability, elec- 
tronic structures, and magnetic properties. We find four 
half-metallic ferimagnetic materials with negative forma- 
tion energies. For the best case of Bi2FeMoC>6, the half- 
metallic gap and Curie temperature reach to 0.71 eV and 
650 K, respectively. This means that high spin polariza- 
tion could be realized at high temperatures well above 
room temperature. More detailed results will be pre- 
sented in the following. 



mented in package VASP— . We use generalized gradient 
approximation (GGA)i^ for the exchange-correlation po- 
tential. In addition to usual valence states, the semicore 
d states are considered for Bi and the semicore p states 
for Cr, Mn, Fe, Mo, and Os. Scalar approximation is used 
for relativistic effect^, and the spin-orbit coupling is ne- 
glected because it has little effect on our main conclusion 
(to be detailed in the following). We use Monkhorst- 
Pack method to generates the K-point meshi£, choosing 
6x6x6 (6x6x4) for structure optimizations and total 
energy calculations of 10-atom (20-atom) unit cells and 
12x12x12 for electronic structure calculations. The cut- 
off energy is set to 500 eV and the criteria for convergence 
is 10~ 6 cV for electronic steps and 0.005 eV/A on atoms 
for ionic steps. 

Metropolis algorithm and variants are used for our 
Monte Carlo simulation s 17 ' 18 . Phase transition tempera- 
tures T c are determined through investigating the average 
magnetization, magnetic susceptibility, and fourth-order 
Binder's cumulant as functions of temperature^. Several 
three-dimensional lattices of upto 30x30x30 magnetic 
unit cells with periodic boundary condition are used in 
these calculations. The first 90,000 Monte Carlo steps 
(MCS) of total 150,000 MCS are used for the thermal 
equilibrium, and the remaining 60,000 MCS are used to 
calculate the average magnetization for a given temper- 
ature. The Curie temperature T c value is determined 
through investigating the magnetization as a function of 
temperature. 



II. COMPUTATIONAL DETAILS 

We use the pseudo-potential and plane wave methods 
within the density functional theory (DFT)iS, as imple- 
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III. MAIN CALCULATED RESULTS AND ANALYSIS 

Comparing Sr 2 FeMoO,A Bi 2 FeCrOe££i&2° and oth- 
ers similar— ~—, we consider double perovskite structure 
of formula Bi2ABC>6, taking some 3d transition- metal 
elements for A and some 4d/5d for B. We optimize 
fully their crystal structures in terms of the unit cell 
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of the 10 atoms. The optimized Bi 2 AB06 has space 
group Rc (#146). This crystal structure, similar to R3c 
(#161), is distorted from cubic double perovskite struc- 
ture, has rhombohedral symmetry, and includes 10 in- 
ternal parameter a 5 i 10 ' 11 . We shall present four of the 
Bi 2 AB0 6 compounds, for AB = FeMo, MnMo, MnOs, 
and CrOs, because they have negative formation energies 
so that their experimental realization should be prob- 
able. Our optimized structural parameters of the four 
Bi 2 ABOg compounds are summarized in Table I. The 
total magnetic moments per formula unit and the partial 
moments in the spheres of the magnetic A and B atoms 
are summarized in Table II. The magnetic moments in 
the spheres of other atoms are much smaller. The to- 
tal moments are integers in unit of Bohr magneton //g, 
showing a feature of half-metallicity^. The magnetic mo- 
ment at the A atom is antiparallel to that at the B atom, 
which means that the magnetic order in these compounds 
is ferrimagnetic. 



TABLE I. Optimized structural parameters of double per- 
ovskite Bi2AB06 with the Rc (#146) crystal structure for 
AB = FeMo, MnMo, MnOs, and CrOs. 
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AB 


FeMo 


MnMo 


MnOs 


CrOs 


a (A) 


5.725 


5.779 


5.761 


5.771 


c(A) 


14.054 


14.093 


13.607 


12.816 


a(°) 


59.91 


60.20 


61.61 


64.37 


Bi 2i 


0.9881 


0.9859 


0.9917 


0.0294 


Bi z[ 


0.4837 


0.4851 


0.4923 


0.5071 


A 22 


0.2589 


0.2572 


0.2593 


0.2682 


B z' 2 


0.7675 


0.7657 


0.7647 


0.2681 


O x 3 


0.5504 


0.5612 


0.5534 


0.4811 


Ox' 3 


0.0506 


0.0471 


0.0520 


0.0558 


O j/3 


0.9357 


0.9220 


0.9273 


0.9353 


oy'z 


0.4350 


0.4397 


0.4348 


0.4287 


O 2 3 


0.1022 


0.1036 


0.0981 


0.1075 


O 2 3 


0.6099 


0.6145 


0.6085 


0.6009 
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FIG. 1. (color online) Spin-resolved density of states (DOS, 
in state/eV per formula unit) of double perovskite Bi2ABOe 
for AB=FeMo (a) and AB=MnMo (b). The solid line is total 
DOS, and short-dashed, dot-dashed, and dotted lines refer to 
partial DOS projected in the atomic spheres of Bi, A, B, and 
O, respectively. The upper part in each panel is majority-spin 
DOS result, and the lower the minority-spin one. 



In Fig. 1 we present spin-resolved density of states 
(DOS, in states/eV per formula unit) between -7.7 
and 3 eV of the double perovskite Bi 2 FeMoOg and 
Bi2MnMoOg. The total DOS in majority-spin channel 
is zero at the Fermi level in both of the cases. This in- 
dicates that the two double perovskite compounds are 
half-metallic, in agreement with the integral magnetic 
moments in unit of /ib- In Fig. 2 we present spin- 
resolved density of states between -8 and 3 eV of the 
double perovskite Bi 2 MnOs06 and Bi 2 CrOs06. They are 
both half-metallic, too, but it is in minority-spin channel 
that the total DOS at the Fermi level is equivalent to zero 
in these two cases. The filled electronic states near the 
Fermi level originate mainly from the B atom (Mo or Os). 
We can use half-metallic gap E g as the key parameter to 
describe the half-metallic property^i^ir— . The E g values 
of the four compounds, from 0.25 to 0.71 eV, are summa- 
rized in Table II. For the Bi 2 FeMo06, E g is equivalent to 
0.71 eV, which implies that high spin polarization could 
be robust even after the spin-orbit coupling is taken into 



We investigate their formation energies to determine 
the stability of these materials towards experimental re- 
alization. For achieving reasonable reliability, we choose 
stable and reachable compounds as our references, and 
try to use those reference compounds whose valence 
states concerned are close to those of our compounds. 
Therefore, we use Bi 2 3 , Cr 3 4 , Mn 3 4 , Fe 3 4 , Mo0 2 , 
and Os0 2 for our reference compounds for calculating 
the formation energies. The formation energy is defined 
as 

E f = £(Bi 2 AB0 6 ) - E Ief) (1) 

where E(X) is the total energy of X, and E le f is defined as 
£(Bi 2 3 ) + ±£(A 3 4 ) + 2£(B0 2 ) - lE{0 2 ). This crite- 
ria is much more severe than merely using AO compound 
or bulk A materials because O atom in the gas state has 
higher energy than in compounds such as Fe 3 4 . This 
should be more precise because the bonds in our mate- 
rials are almost formed between metal atom and O, not 
between metal atoms. The formation energies for the 
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FIG. 2. (color online) Spin-resolved density of states (DOS, 
in state/eV per formula unit) of double perovskite B12ABO6 
for AB=MnOs (a) and AB=CrOs (b). The solid line is total 
DOS, and short-dashed, dot-dashed, and dotted lines refer to 
partial DOS projected in the atomic spheres of Bi, A, B, and 
O, respectively. The upper part in each panel is majority-spin 
DOS result, and the lower the minority-spin one. 

four compounds are summarized in Table II. The nega- 
tive values means that they should probably be realized. 



cells. Rigorously speaking, there are some induced spin 
density in the spheres of the Bi and O atoms, less than 
0.05^b- Because they are very small compared to those 
in the spheres of the magnetic atoms, we shall consider 
only the magnetic atoms in the following. Actually, A 
and B atoms form a lattice of magnetic unit cells (cu- 
bic unit cells of a NaCl crystal structure) ^ 10 ' 11 . In these 
calculations, we fix the structures and change the mag- 
netic orders of A and B atoms. In order to make the 
electronic steps converge for a magnetic order, the linear 
mixing parameter should be decreased to an small value, 
0.1 or smaller. Through comparing the total energies, we 
obtain the spin exchange interaction constants: Jab for 
the nearest A and B pair, Jaa and Jbb for the A-A and 
B-B next nearest pairs. Jab is dominant over the others. 
The resultant spin Hamiltonian reads: 

H = ^ JijSi ■ Sj (2) 

to) 

where Si is spin operator at site i (in both of the A and 
B sublattices), the summation is over spin pairs, and the 
spin interaction constant Jy is limited to the nearest and 
the next nearest neighboring spins. 
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TABLE II. Calculated values of formation energy (Ef in eV), 
magnetic moment of atom A (Ma in ^b), magnetic moment 
of atom B {Mb in /is), total magnetic moment (M in \ib 
per formula unit), half- metallic gap (E g in eV), and Curie 
temperature (T c in K) of double perovskite Bi2AB06 for AB 
= FeMo, MnMo, MnOs, and CrOs. 



AB 


FeMo 


MnMo 


MnOs 


CrOs 


M A 


3.638 


4.279 


4.066 


2.636 


Mb 


-1.755 


-1.391 


-1.041 


-0.613 


M 


2.000 


3.000 


3.000 


2.000 


E a 


0.71 


0.47 


0.46 


0.25 


Ef 


-0.41 


-0.26 


-0.42 


-0.29 


T c 


650 


255 


174 


201 



In order to estimate the Curie temperatures (T c ) of 
the materials, we calculate the spin exchange interac- 
tions between the nearest and next nearest neighboring 
magnetic atoms (A and B) in terms of the 20- atom unit 



FIG. 3. (color online) Average normalized magnetizations as 
functions of temperature for double perovskite Bi2FeMoOe 
for four different L values. Monte Carlo simulations are done 
with L x L x L magnetic unit cells. 

We carry out Monte Carlo simulations to estimate the 
T c of the material o 17 ' 18 . It is well known that Curie 
temperature will be a little underestimated if classical 
approximation to the Heisenberg model (2) is used in 
the Monte Carlo simulation, but it must be much over- 
estimated if the model (2) is reduced to Ising model. 
For comparison, we do our Monte carlo simulations with 
both of the approximate models. We present the aver- 
age normalized magnetization of Bi2FeMo06, from clas- 
sical Heisenberg model, as a representative in Fig. 3. 
The T c value can be estimated to be 650 K. The others 
can be done in the same way. The calculated T c val- 
ues are summarized in Table II. In contrast, the Ising 
model results are 1010, 396, 264, and 270 K for AB = 
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FeMo, MnMo, MnOs, and CrOs, respectively. There- 
fore, the Curie temperatures for the four half-metallic 
ferrimagnets are, at least, 650, 255, 174, and 201 K for 
Bi 2 FeMo0 6 , Bi 2 MnMo0 6 , Bi 2 MnOs0 6 , and Bi 2 CrOs0 6 , 
respectively. High Curie temperature well above room 
temperature could be realized in Bi 2 FeMoC>6. 

IV. DISCUSSIONS AND CONCLUSION 

Our calculated results shows that the spin exchange 
interaction between the nearest A and B atoms is posi- 
tive, and the A-A and B-B interactions are either weak 
or negative depending on specific A and B atom o 10 ' 11 . 
In the case of Bi 2 FeMo06, our calculations show that 
the nearest A-B spin exchange energy is 39.2 meV, and 
the nearest A-A and B-B spin exchange energies are 0.13 
and -0.71 meV, respectively. The main spin interaction 
is intermediated by the O atom in between the magnetic 
A and B atoms, with the A-O-B bond angle being al- 
most 180°, and therefore, it is an antiferromagnetic su- 
perexchange. The A atom contributes a different mag- 
netic moment from the B atom so that ferrimagnetism is 
formed in these double perovskite compounds. Possible 
overlapping of the nearest O wave functions should play 
some roles in these compounds, but the main mechanism 
for the ferrimagnetism must be the antiferromagnetic su- 
perexchange between the nearest A and B atoms. 

In summary, our first-principles calculations show that 
four double perovskite oxides, Bi 2 ABC>6 (AB = FeMo, 
MnMo, MnOs, and CrOs), have negative formation en- 
ergy, from -0.42 to -0.26 eV per formula unit. In the case 
of Bi 2 FeMo06, our calculated results uncover that its 
half-metallic gap and Curie temperature reach to 0.71 eV 
and 650 K, respectively. These indicates that they could 
probably be realized and high spin polarization could be 
achieved at high temperature. We believe that at least 
some of them could be synthesized soon and would prove 
useful for spintronic applications. 
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